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Abstract 

We study how to construct Dirac fermion defined on the honeycomb lattice in position 
space. Starting from the nearest neighbor interaction in tight binding model, we show 
that the Hamiltonian is constructed by kinetic term and second derivative term of three 
flavor Dirac fermions in which one flavor has a mass of cutoff order and the other flavors 
are massless. In this formulation the structure of the Dirac point is simplified so that its 
uniqueness can be easily shown even if we consider the next-nearest neighbor interaction. 
We also explicitly show that there exists an exact chiral symmetry at finite lattice spacing, 
which protects the masslessness of the Dirac fermion, and discuss the analogy with the 
staggered fermion formulation. 
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1 Introduction 



Graphene forms from a layer of carbon atoms with hexagonal tiling pQ and it is much 
discussed in not only condensed matter physics but also high energy physics for its remarkable 
features [2]. One of the most important features of graphene is that the quasiparticles behave 
like massless Dirac fermion with effective speed of light near c/300 [3]. An explanation to 
the question why massless Dirac fermion emerges in non-ralativistic many body system was 
primarily given by Semenoff 0]. Ref. [1] has shown that the low energy excitations around 
two independent Dirac points on the fermi-surface can be described by two relativistic Weyl 
fermions having opposite chiralities, which are also regarded as massless Dirac fermion. 
However, since this construction is based on the momentum space representation, it is not 
clear how one can introduce the effect of gauge interaction as local interactions. 

In this paper we show the construction of Dirac fermion from honeycomb lattice in posi- 
tion space, which allows one to introduce local gauge interactions and perform the dynamical 
calculation of physical observables. Furthermore in this formulation the structure of flavor 
symmetry becomes manifest, and we can easily identify the quantum number of low energy 
excitations. We also expect that the formulation enables to advance a study of the dynami- 
cal nature of graphene such as the gap generation, which is associated with the spontaneous 
breakdown of chiral symmetry. 

This paper is organized as follows. In section [2j we briefly review the construction of 
Dirac fermion from honeycomb lattice in momentum space. In section [3j after introducing 
the formulation defined in position space, we discuss uniqueness of Dirac point and physical 
modes, and then we show that two massless Dirac fermion appear at low energy region in 
section 4. In section 5, we discuss a derivation of exact chiral symmetry in our formulation, 
and the last section is devoted to the summary and discussion. 

2 The conventional derivation from honeycomb lattice 

We first review the conventional derivation of Dirac fermion formulation from tight binding 
model of honeycomb lattice jH |8] . Let us start from the tight binding Hamiltonian 

r&A i=l,2,3 cr 

[ E EE 4(r>.(r + s ;o + E E E bt ^+ fyw] - 0) 

r£Aj=l cr r£B j=l cr 

where the first line is a nearest neighbor hopping term and the second line is a next nearest 
neighbor hopping term, and t,t' are hopping amplitudes, a), a and b\b are the fermionic 
creation and annihilation operators of electrons on two triangular sublattices A and B in 
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Figure 1: Honeycomb lattice is constituted of two triangular sub-lattices A and B, which 
are colored with red and blue respectively. 



Figure |l| Si(i = 1,2,3) and s'-{j = 1, • • • ,6) denote the relative position vectors for three 
nearest neighbors, and for the six next-to-nearest neighbors respectively, and Si{i = 1,2,3) 
are given explicitly by 

si = a ( 1, ) , s 2 = a ( -1/2, ^3/2 ) , s 3 = a ( -1/2, -y/3/2 ) , (2) 

where ao is a lattice spacing. We note that in graphene system above parameters are given 
as t = 2.8 eV and if = 0.1 eV [9], and ao = 1.42 A [2]. In the following discussion we take 
t' = and also suppress spin index a. In order to consider Dirac point on fermi-surface, the 
Fourier transformation 



a(r) 



d 2 k 



ik-r 



a(k), b(f) 



d 2 k 



ik-r 



b(k) 



(3) 



(2vr) 2 v " VJ J (2vr) 2 
for the fermionic creation and annihilation operator is usually employed, and thus the Hamil- 
tonian in momentum space is given by 



n 



d 2 k ( a(k) 
W? V b{k) 



D(k) 
D*(k) 



o(Jfe) 
b(k) 



with 



D(k) = t ^ e 



i=l,2,3 

The energy eigenvalue of the above Hamiltonian is represented as 

E(k) =±t\ eilsl 

j=l,2,3 



(4) 



(5) 



(6) 
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In the half filling electron system, the negative and positive eigenvalues corresponding to the 
valence band and conduction band respectively appear, and it turns out that there exist two 
independent Dirac points K± on the fermi-surface, which satisfy E{K±) = 0. 

In order to derive the low energy effective Hamiltonian, we expand D(k) around the Dirac 
points to first order of momentum p. Regarding K± and A, B site as spin degrees of freedom 
(DOF), and defining 4-components spinor as 

m = {a(K++p),b(K++p),b(K- +P),a(K- +p)f, (7) 
the Hamiltonian in Eq.Q reads 

n w w E / (02 f (p) Moh\ • (8) 

Here v = 3aot/2 is interpreted as a fermi velocity of quasiparticles, and the matrices 70, 71, 72 
satisfy Clifford algebra {7^,7,} = g pu • l4x4> where g pu is a metric in 2 + 1 dimensional 
space-time. Further, introducing additional matrix 73, which anti-commutes with the above 
70 j 7i) 72; we can define 75 = ^0717273, which is a generator of chiral symmetry forbidding 
a (parity-invariant) mass term m^jo^. 

In the next section we will show an alternative derivation of Dirac fermion based on 
position space, which has manifest locality and chiral symmetry. 



3 Formulation in position space on honycomb lattice 
3.1 Tight binding model on the real space lattice 

In this section, we explain our formulation and address Dirac point and physical modes. First 
we consider the new labeling of DOF of the fermionic creation and annihilation operator as 
shown in Figure [2J In this labeling we define Ap and Bp (p = 0,1,2) as the new DOF 
having the operators on the site of honycomb lattice and the argument of operator location 
is defined as the position of center of hexagonal unit cell with position vector e p (p = 0, 1, 2). 
Thus in our formulation the tight binding Hamiltonian is expressed as 

n = ST^sr^( XA P (x) \ / t'U(x,y) pp/ t<f>(x,y) pp < \ / XA P '{y) \ (9) 

hb V ) V J V xbAv) )' 
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Figure 2: Honeycomb lattice 



where <&(x,y) and IL(x,y) are given by 

f T 1 1 \ 



y) 



U(x , y) 



\ 



1 Ti 1 
1 1 T 2 I _ 

' X 



1 + 4 + :/, 



■ u 



1 + T q + t\ 1 + To + T\ \ 



\ 1 + T\ + T 2 1 + T[ + T 2 



1 + T x + 




(10) 



(11) 



with 



(Tp)x,y 3x,y+e p i (l)z,j/ ^x,y- 



(12) 



Here x, y denote lattice sites on the fundamental lattice, which is composed of hexagonal 
unit cells bounded by red circles in Figure [2j and basic vectors e p (p = 0, 1, 2) are given by 



e = a 1, 



ei = a 



1/2, ^3/2 , e 2 = a( -1/2, -^3/2 



(13) 



where a is a lattice spacing as a distance between hexagonal unit cells and these vectors 
satisfy eb + ei + e 2 = . xl p (^)i Xip{%) are the new defintion of creation and annihilation 



d Here we note that there is following relation between e p (p = 0, 1, 2) and Si(i = 1, 2, 3): 

e = 3s*i, ei = 3s 2 , e 2 = 3s*3- 



(14) 
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operators on the triangular sublattice /(= A, B) of honeycomb lattice. And we note 

n(f, y) = J2 ®& y) ~ 3 • ^xs (is) 

z 

holds, which will be useful later. 



3.2 The nearest neighbor interaction 

In this section we assume t' = as the leading tight binding approximation, and we discuss 
the general representation of Eq.Q. 

Let us introdce the first and second derivatives on the lattice, 

VpX(x) = ^(x{x + e p ) -x(x-e p )), (16) 

ApX(x) = X(x + e p ) + x{x - e p ) - 2x(x). (17) 

The leading order of Hamiltonian Q is written as the following tensor product form; 



(n ® M)x(x) - i ^(r 2 T p )(V pX (x)) + - £(n ® r p )(A pX (x))j (18) 



with 



and 



X(«) = (XA0(^),XAl(^),XA2(^),XB0(^),XBl(^),XB2(^)) rj 
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1 
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1 
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/ 
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1 






r = 















, r 1 = 







1 







, r 2 = 

















V 
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V 











J 
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J 




V 
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/ 



(19) 



(20) 



It is clear that the first derivative term of Eq.(18) corresponds to kinetic term, and the sec- 



ond derivative term corresponds to flavor dependent Wilson-like term. 



Next we discuss a constant mode which is the eigenmode of the first term of Eq.(18) 

M = tJ2xHx)(n^M)x(x) 



t 



tM, 



W( XAp ^ \ ( U tMpp ' \ ( Xa p'^ \ 
sh\ W*) J V tM pP > ) { XBp ,(x) ) ■ 

Using the diagonalization of mass matrix tM, the above equation is described as 



i>A P {x) 



M diag 



x pj \ ^B P {x) 

f St \ 







(M di ^) pp , 
(M di ^) pp , 



ipAp'(x) 
*pBp>(x) 



(21) 



(22) 



(23) 



with unitary transformation, 

Xl P {x) 



7s S ' 

V P'=0,l,2 



,i2irpp' /3 



(24) 



and thus we interpret from the above diagonalization that the constant mode can be decom- 
posed into two zero modes and one massive mode. 

In order to prove the uniqueness of the Dirac point, we consider an eigenvalue equation 
of this Hamiltonian in momentum space. Using following Fourier representation of xip(cc), 

d 2 k 



Xi P (x) 



Ak-x~, 



(2*) s 



Xip(k), 



(25) 



Hamiltonian ( 18 ) in momentum space is given by 
d 2 k 



(27T)' 



E 

p,p' 



XAp(k) 
XB P (k) 



t 











XA P '(k) 
XB P '(k) 



(26) 



where 



$(fc) 



( bo 1 1 \ 

1 bi 1 
V 1 1 b 2 j 



(27) 



with bp = exp [— • e p ] (p = 0, 1, 2). The eigenvalue equation can read into 



"A 2 • 1-3x3 



(28) 



where 



*(jb)*(fc) 



(29) 



/ 



3 b + b\ + 1 6 + ^ + 1 \ 

+ + l 3 6i + 6^ + l 

\b* + b 2 + l b\ + b 2 + l 3 

Assuming the electron system is half-filling in which the Dirac point on fermi-surface should 
appear at A = 0, we can derive the condition, which Dirac point fulfills, by the substitution 
of A = into the above eigenvalue equation; we solve an equation |$(fc)$t(fc)| = Using 
606162 = 1 (since e*o + e*i + e*2 = 0), we find a simple condition 



Zk 



In general the eigenvalue equation of Eq. ( 28 1 can be represented as 

F{\ 2 ) = 0, F(x) =x i - 9x 2 - 3(z fe + z* k - 6)x - \z k - 3j ; 



(31) 



(30) 



where Zk will be defined in Eq. 1 32 1 . 
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where 



Zk = b + h + b 2 = e 



+ e 



-ik-e-i 



+ e 



-ik-&2 



(32) 



and from this condition we have 



-ik-eg 



-ik-e-i 



-ik-e2 



1, 



(33) 



which means k = n\U\ + n2'U2^](ni, n 2 : integer), therefore we prove that Dirac point appears 
only at k = within the BZ. It turns out that the constant mode provides unique Dirac 
point at k = in this formulation within the nearest neighbor hopping approximation in 
contrast to the traditional formulation in which there are two Dirac points at the edge of BZ 
(see the section [2|. 



3.3 The next-to-nearest neighbor interaction 

Here we will confirm that the Dirac point explained in the previous subsection does not 
disappear if we add the next-to-nearest neighbor hopping term. Turning on t' 7^ 0, the tight 
binding Hamiltonian (|9|) is written as 



U 



EE 



x,y p,p 



XAp(x) 
XBp(x) 



- 3< ee( 



XAp{X) 

XBp(x) 



1 3 X3 
1 3 X3 



XAp'(y) 

y XBp'iv) 

XAp'iy) 
XB P '(y) . 



(34) 



where we substituted Eq.(15). The second line merely shifts the origin of energy, and does 



not affect the dynamics at all. Therefore we neglect the constant term in the following 
discussion. 

When we naively consider constant mode as in the previous subsection, the mass matrix 
of this Hamiltonian becomes 



3? 



1 1 1 



\ 



V 



V 



1 1 1 
1 1 1 

(l 1 l\ 

1 1 1 
1 1 1 



1 1 1 



/ 



3t' 



1 1 1 

1 1 1 

(l 1 l\ 

1 1 1 

1 1 1 



(35) 



/ 



f Note that {ui}i=i,2 is the basis vector of reciprocal lattice space, which is defined as e*i ■ itj — 2n8ij 

(*,i = i,2). 
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then transforming the same basis as ipi p (x) in Eq.(25), this matrix is given as 



V 



{ 9t \ 





/ 3t \ 







( St \ \ 





( 9i' \ 







/ 



/ 



(36) 



and thus two zero modes and one massive mode also appear. One may think this implies that 
Dirac point does not change even if the next-to-nearest neighbor hopping term is included. 
However, it is not so trivial whether the fermi-surface remains in vanishing energy level 
holding on the half-filling configuration. 

In order to investigate in detail, we discuss the energy spectrum of the modified Hamil- 
tonian in momentum space, which is expressed as 

d 2 k 



(2vr)< 



xHk)[K(k)]m 



with 



K(k) 



t'$(k)&(k) t$(k) 
t&{k) t'&(k)$(k) 



We diagonalize $(fc) by the following bi- unitary transformation 

U A ${k)U B = diag( 



(37) 



(38) 



(39) 



where Ua,Ub are some unitary matrices and <f>i (i = 0,1,2) are real eigenvalues of <&(&). 
Performing unitary transformation U = Ua © Ub, K{k) is transformed into 



K{k) 











t(f>o 











t'rt 








% 



















t<h 


t<f>o 



















t(f>l 



















t<h 











(40) 



where $ {% — 0,1,2) are the eigenvalues of $(/s)$t (k), and satisfy the eigenvalue equation 

F{tf) = 0, (41) 
where F{x) is defined in Eq.(|30j). Here, the eigenvalue equation of K{k) is expanded into 



n k*' 



t 



t 



2 j.21 



0, 



(42) 



i=0,l,2 
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where A' is an eigenvalue of K(k), and we obtain 



£4(&) = t>?±t|&| 

= t'4>lf±{\<j>i/<j> c \) (43) 



with (f) c = t/t' and f±(x) = x 2 ± x. In Eq.(43), taking t' = 0, the fermi-surface appears at 
vanishing energy level in the half-filled electron system since the number of positive eigen- 
values E + ((f>i) = t\4>i\ and the number of negative eigenvalues E_(4>i) = — t\(pi\ are the same. 
However, taking account into the effect of the next-to-nearest neighbor hopping term, t' ^ 0, 
the situation turns less simple. 

We plot E'±((j)i) jt '4> 2 C against \4>i/4> c \ in Figure [3j It is obvious that the number of eigen- 
values E' + ((f>i) and that of E'_{(j)i) are also the same. In Figure |3j we see that the negative 
eigenvalues E'_{4>i) remain in negative values unless \4>i\ exceeds \(f> c \, while the eigenvalues 
E' + (4>i) stay in positive values at arbitrary \<fii\. Thus if \<fii\ does not exceed the threshold, 
the fermi surface remains in vanishing energy since the number of the positive eigenvalues 
equals that of the negative eigenvalues. On the other hand, if \<f>i\ exceeds the threshold, the 
fermi surface stays no longer in the same energy level. 

From the eigenvalue equation of <&(k)& (k) in Eq.(41), the range of \(f>i\ is shown to be 
between and 3. Therefore, when 3 < \(f> c \, the Dirac point is stable in adding the next- 
to-nearest neighbor hopping term. Especially in graphene system, the fermi surface is not 
changed since \cp c \ ~ 28. 



4 Close to the continuum limit 

In this section, we show the derivation of Dirac fermion formulation in position space close 
to the continuum limit. Here we consider the low-energy expansion of Hamiltonian around 
Dirac point and investigate "chiral" symmetry which forbids mass gap in low energy region. 

In momentum space, the tight binding Hamiltonian in terms of the mass eigenstate is 
given as 




where ipi a {k) is a Fourier representation of the mass eigenstate of ipi a (x) and $'(A;) denotes 
<£(£;) in basis of these mass eigenstates, which are defined as 

^Ia(k)=J2^'^Ia(x), (45) 

x 
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-1.0 - 



Figure 3: x axis denotes \cj)i/<p c \, and y axis denotes energy eigenvalues divided by t'4> 2 . The 
dashed lines colored with red and blue are E + ((pi)/t'(p 2 = \(pi/(j) c \ and E-(c/>i) /t' '4> 2 C = —\4>i/4> c \ 
respectively. While the solid lines colored with red and blue are E' + ((pi)/t'(j) 2 = f + (\4>i/(j) c \) 
and = f-ilb/M)- 



^ 6 + 61 + 6 2 + 6 6 + uj 2 bi + ub 2 b + ujbi + u?b 2 ^ 

bo + u)b\ + uj 2 b2 6o + 6i + &2 — 3 60 + uj 2 b\ + ub2 
y 6 + uj 2 bi + ujb 2 6 + cj6i + uj 2 b 2 6 + 61 + 6 2 - 3 J 



(46) 



Expanding with respect to momentum around Dirac point to the second order in k as 

-i i \ 

1 -i 
-i i 
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/ 
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1 
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V 



I 



( 



a 

T 



+ 0(A; 3 ), 



(47) 



fc 2 + A| (^ + ^2)2/2 (fc! 

(ki-ik 2 ) 2 /2 kf + kj (k 1 + ik 2 ) 2 /2 
\(h + ik 2 ) 2 /2 {kx-ik 2 ) 2 /2 k 2 + k 2 

where ki(i = 1,2) are components of momentum in Cartesian coordinates, and integrating 
out the massive mode ipio(k) by its equation of motion 



/ 



*'(A:)oo \ ( ^Ao(k) 

&(k)00 J \ ^Bo(fc) 



E 

a=l,2 



*'(A:)o« \ / i>Aa(k) 



(48) 
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the tight binding model is reduced to the following form to the second order in k: 

n cS = J -^^\k)[H x {k) + H 2 {k)]^{k), (49) 

$(k)=($Ai$) Mk) $ m {%) 4>B2(k)) T , (50) 

Hi{k) = v^ki (r 2 ® cri) + k 2 (r 2 (g) ct 2 )}, (51) 



# 2 (A0 = ^| 2 ( fc f + Wi ® 12x2) + " Wi ® *i) - (*ifca)(n ® <r 2 ) j (52) 

with v = at/2, where the former of tensor structure acts on sub-lattice space I = A,B 
and the latter of tensor structure acts on flavor space a = 1 , 2 in the above tensor product 



representations. Here the first term of Eq.(49) is Dirac Hamiltonian in momentum space 



while the second term of Eq.(49) is a higher derivative term suppressed by the mass of ipio, 
which is negligible in the low energy region. 

Close to the continuum limit, we can neglect higher derivative terms H 2 (k) in Hamiltonian 



(49), and the effective Hamiltonian has following 4 global symmetries in this approximation: 

1-2x2 ® 1-2x2, Tl®a 3 , T 2 (g>l 2x2 , T 3 <g> (J 3 . (53) 

Considering a Parity invariant mass term, which is invariant under the Parity transformation 
(which is exchange symmetry of A o B and x — > — x but p — > p in Figure [2]) , 

mi)\k){n® 12x2)$$), (54) 

two global symmteries, l 2x2 ® l 2 x2, t\ ® 03, are preserved, whereas t 2 ® l 2x2 , 73 ® 03 are 
broken. We call the symmetry under global rotation with generators t 2 (g> l 2x2 , T3 <g) 03 as 



"naive chiral symmetry", which forbid the Parity invariant mass term Eq.(54), up to the 
first order of A; (in Table [lj. 



However the higher derivative term, which is the second term in Eq.(49), violates "naive 
chiral symmetry" like Wilson fermion formulation. This fact seems that the mass term may 
be easily induced through quantum corrections when we introduce interactions between elec- 
trons. In general, symmetries at finite lattice spacing may be deformed by lattice artifact. 
We consider a possibility of modified chiral symmetry which depends on the momentum and 
determine its form exactly in next section. We will find that there exists a chiral symmetry 
which is reduced to the above chiral symmetry T3 (g) 03 in the continuum limit. 

Note that in Legendre formulation Lagarangian is descirbed as 

£ = itp(t,x) <9o7o - v ^2 "Yidi VK^) (55) 

i=l,2 
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where ip = V^7o and gamma matrices 70, 71, 72 are, 

7o = ( ° 1 ) , 7i = ( mi J , 72 = ( ^ J (56) 

(in details see also appendix) . Apparently, these gamma matrices 70 , 71 , 72 satisfy Clifford 
algebra {7^,7^} = guv, where is a metric in 2+1 dimensional space-time. These gamma 
matrices are connected 7^ in section 2 by an appropriate unitary transformation. 

Table 1: Symmety in Effective Theory to the mass term 



Global symmetry 


preserved 


broken 


12x2 ® 12x2 


n ® fj 3 


T2 ® 12x2 


T3 ® 0-3 



5 Exact chiral symmetry 

In this section we consider the exact chiral symmetry in position formulation of graphene 
model under leading order of tight binding approximation. We consider a modified chiral 
symmetry at finite lattice spacing involving the first and second derivatives defined in Eq. ^ . 



Here we take following ansatz for exact chiral symmetry of Hamitonian (18) in position 
space, 



5x(x) 



[t 3 ® X) X (x) + \ Y,(T3 ® Y p )(A pX (x) + 2xOH)) + 7 I> ® ^)(V P x(^))l , (57) 



where X, Y p , and Z p are 3x3 hermite matrices. Based on this ansatz, we look for the solution 
for the symmetry equation 5% = 0, and finally we determine its concrete form uniquely |^} 
An explicit calculation of X, Y p , and Z p is shown in appendix [Bj and here we show the result 
as follows: 



5 X (x) = ieh(r 3 ® X) X (x) + \ ]T(T3 ® Y p )(A pX (x)) + - ® Z p )(V pX (^)) 



(58) 



Generally, there are two possibilities for Ts in Eq.(69l, which is consistent with the results in section 



4. One is F 5 = t$ <g> Xi + O(a), which is supposed above derivation, and the other is f 5 = T2 ® X2 + O(a), 
which might be consistent with naive chiral symmetry T2 ® 12x2 in effective theory. Thus, we consider both 
cases. However the latter is found not to satisfy the symmetry equation at the second order in momentum 
expansion. 
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(59) 



(60) 



(61) 



If we take a continuum limit a — > 0, above chiral transformation becomes 5x{x) = 3i8Xx(x). 
Here, taking mass basis, X may be transformed to following form, except for an overall factor, 



(62) 



Thus, the exact chiral symmetry corresponds to chiral symmetry in the continuum T3 (8) o 3 
in effective theory, whch is consistent with the results in section 4. 
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6 Summary and discussion 

In this paper we study how to construct Dirac fermion from honeycomb lattice in position 
space. In our formulation we make new labeling of fermion field in which the fundamental 
lattice is composed of the centers of hexagonal unit cells. The six sites in each unit hexagonal 
cell is reinterpreted as spin-flavor degrees of freedom. 

Starting from the tight binding model, we find that the Hamiltonian is constructed by 
kinetic term and second derivative term of three flavor Dirac fermion in which one flavor has 
a mass of cutoff order and the other flavors are massless. In our formulation the structure of 
the Dirac point is simplified so that its uniqueness can be easily shown. We also explicitly 
show that there exist an exact chiral symmetry at finite lattice spacing, which protects the 
masslessness of the Dirac fermion, under the nearest neighbor interaction. 

Prom the point of view in lattice gauge theory, the graphene model introduced in ref . [4J 
or in this paper seems analogous to the staggered fermion formulation [5|. In the staggered 
fermion formulation, there exists two representations of Dirac fermion with grassmann vari- 
ables spread over hypercubic lattice. One is defined in momentum space [6] and the other 
is defined in position space [7]. In the former formulation the momentum space is divided 
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into 2 d regions (where the size of the Brillouin Zone (BZ) becomes half in each direction) 
and the staggered fermion field defined in these regions can be regarded as 2 d / 2 flavor Dirac 
fields. While, in the latter case, the staggered fermion defined in hypercubes which has 2 d 
corners is introduced and Dirac field with 2 d spin-flavor is defined by multi-product of Dirac 
gamma matrices on grassmann variables. Comparing to the Graphene model suggested by 
Semenoff, in which Dirac fermion is defined as two excitations on different regions in BZ, it 
is very analogous to the staggered fermion in momentum space formulation [Bj. On the other 
hand, our formulation is analogous to the position space formulation for staggered fermion 
found in ref.jT]. 

The formulation in position space easily extends toward the gauge interacting system and 
also provides the complementary ingredients for understanding of the connection to QED 
with 2+1 dimensional fermion brane |1CH lllj . Monte-Carlo simulation of non-relativistic 
action \\.2\ 1X3] and honeycomb lattice simulation |144 115] . Furthermore, since our formula- 
tion has manifest structure of flavor symmetry, the analytical study of dynamical nature of 
graphene such as the spontaneous symmetry breaking in this formulation will be comparable 
with lattice simulations. 
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A Model in Lagrange formulation 

If we add a mass term to the effective theory, Hamiltonian is written as follows; 

/d 2 k 
7 —^^(k)[a 1 ki + a 2 k 2 + mP]tp(k), (63) 
( 27r ) 



where 



ai = ° V a 2 = ° "** ]. (64) 

V Utx ) I ia 2 



j3 is a hermitian matrix and we may take following choices; 

(65) 



1 \ I 

1 } ' I -1 



Here the first gives parity even mass term and the second gives parity odd mass term. 
However the parity odd mass term may be forbidden by parity, thus we choose the parity 
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even mass term here. Then transforming above Hamiltonian to Lagrangian in real space, we 
obtain following Dirac Lagrangian in configuration space. 

C = iip\t,x) do + v otidi — m(3 ip(t,x) 

i=l,2 

= iip(t, x) (9 7o ~ v ^ lidi ~ m ip(t, x) 

i=l,2 

where ?/> = and 70 = (3, 71 = — /3ai, 72 = — /3a2- Evidently the gamma matrices 
70) 7i)72 satisfy Criford algebra {7^,7^} = g pu , where g pu is a metric in 2+1 dimensional 
space-time. 



(66) 
(67) 



B Explicit calculation of exact chiral symmetry 

In order to determine X, Y p , Z p , we employ momentum representation of x(x), X^{x 
d 2 k 



H 



I 



(27T) 



2 X f (fc) [(71 ® A) + ^ e^(r_ ® T p ) + e~ lkp (r + ® T P )J (68) 
p p 

with T-t = (n ± U2)/2 and A = M — 1, and for chiral transformation <5x(&) = z0rs(A;)x(A;) 
f^fc) is given as 



f 5 (fc) = (r 3 ®I)+^ e*Sp + £ e ~^P' 



with 



Here, imposing [H(k), T^(k)] = 0, we obtain following equations; 

{A,x} + £(r,w, + wJr p ) = o 
p 

{r p ,X} + AWt + W p A = 

ah/ p + wJa + £ owj + w A rv) = 

r P wJ + H/ P r p = 

r p w a + w^r p = o (p^a). 

Solving these equation for X and W p (p = 0, 1, 2), the solutions are found to be 



X 



-i i 

1 -i 
-i i 
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(69) 

(70) 

(71) 

(72) 
(73) 

(74) 
(75) 

(76) 
(77) 
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